Multiscale entanglement in ring polymers under spherical confinement 
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The interplay of geometrical and topological entanglement in semiflexible knotted polymer rings 
confined inside a spherical cavity is investigated using advanced numerical methods. By using 
stringent and robust algorithms for locating knots, we characterize how the knot length Ik depends 
on the ring contour length, L c and the radius of the confining sphere, R c . In the no- and strong- 
confinement cases we observe weak knot localization and complete knot derealization, respectively. 
We show that the complex interplay of Ik, L c and R c that seamlessly bridges these two limits can be 
encompassed by a simple scaling argument based on deflection theory. The same argument is used to 
rationalize the multiscale character of the entanglement that emerges with increasing confinement. 

If we tie a knot in a piece of rope and pull the latter at its extremities, the knotted part will become readily 
distinguishable from the rest of the rope because it localizes. Indeed, for a very long rope only a neglible fraction of 
its contour length will be required to accommodate the knot pQ. In this intuitive example of knot localization the 
topological (global) entanglement embodied by the knot does not interfere with the geometrical (local) entanglement 
of the rest of the rope. This is hardly the case for polymers that circularize in equilibrium. In such rings, knots are 
abundant [5] and the interplay between topological and geometrical entanglement impacts significantly the molecules' 
dynamical, mechanical and metric properties [3HI] . Understanding how and to what extent global and local entangle- 
ment are related is a major open issue in polymer physics [B] with ramifications in key biological contexts[7|, especially 
those related to genome organization in eukaryotes, bacteria and viruses [Br[T2"]. 

A first breakthrough in the problem could be made by establishing what fraction of the polymer ring is occupied 
by the knot(s) and whether this measure is sufficiently robust or depends instead on the geometrical entanglement 
degree. As a prototypical context to examine this problem we consider semi-flexible, self-avoiding rings of cylinders 
with the simplest knotted topology, a 3i (trefoil) knot, and subject to isotropic spatial confinement. In such a system, 
by varying the size of the confining region, the degree of geometrical complexity can be changed and related to the 
equilibrium size of the knotted ring portion. 

For definiteness, the rings properties are set to match dsDNA. Specifically, the cylinder diameter and long axis are 
set equal to d = 2.5nm and b = 10nm, respectively. The latter quantity is ten times smaller than the DNA Kulm 
length (equal to twice the persistence length, l p = 50nm) thus ensuring a fine discretization of the model DNA. The 
system energy includes steric hindrance of non consecutive cylinders plus a bending potential: 

/ N 

i=i 

where ti is the orientation of the axis of the ith cylinder, t^+i = t\ and the temperature T is set to 300K. We 
considered rings of N = 50, ...250 cylinders, corresponding to countour lengths, L c = Nb ranging from 500tito to 2.5 
fj,m. This range allows for probing changes in knot localization going from semiflexible to fully-flexible rings |13j as 
well as examing the effect of the interplay between L c , l p and the radius of the confining sphere, R c . For simplicity 
of notation in the following all lengthscales are expressed in units of b. 

Because compact ring configurations are entropically disfavoured with respect to unconstrained ones, simple stochas- 
tic sampling schemes cannot be effectively used to generate spatially confined rings [6 . An analogous entropic attrition 
works against having a sizeable population of knots of a given type, trefoils in our case, at all levels of confinement [S]. 
To overcome these two difficulties we used a biased Multiple-Markov-chain sampling scheme [12] with 24 Markovian 
replicas. For each replica, ring configurations are evolved using crankshaft and hedgehog Monte Carlo moves Q3]. The 
moves preserve the length of the rings, but not the topology, consistently with ergodicity requirements [6, 15 . In fact, 
even when the configurations before and after the move are self-avoiding, the move itself may involve self-crossings 
of the chain. A newly generated ring, T, is accepted according to the standard Metropolis criterion with canonical 
weight, exp [— U(T)/ KgT], where the potential energy, U is suitably chosen to enhance the populations of compacts 
trefoil rings. Specifically, if T is non self-avoiding or has a too complicated topology (8 or more minimal crossings 
according to the Alexander polynomial [6]) then U is set equal to infinity, otherwise U(T) = E{,(T) +pR c (T). In the 
latter expression, R c is the radius of the sphere containing T (calculated as the largest distance of the cylinder vertices 
from the ring geometric center) and p is a confining field whose value is set differently for the various replicas so to 
sample rings with different degree of confinement. 
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FIG. 1. (a) Typical configurations of trefoil-knotted rings with N = 200 cylinders at different levels of confinement. The 
encapsidated configurations for R c = 6 and R c = 3 are magnified for ease of visualization. The knotted portion of the rings is 
shown in red. (b) Sketch of a trefoil-knotted ring; its knotted portion is highlighted in red in panel (c) and its shortest knotted 
arc is shown in blue in panel (d). (e) Average knot length, (t, as a function of 1/R C for rings of varying number of cylinders, 
N = 50. ..250. (f) Dependence of Ik on L c in the no- and strong-confinement limits; the lines represent the best-fitting power 
laws. The strong-confinement data correspond to the limiting values approached exponentially by the curves in panel (e). The 
exponential fit was limited to the data points for which Ik > 0.6L C . 



For each value of N this procedure was used to generate at least ~ 10 7 configurations with varying compactness 
and topology. For the a posteriori data analysis we picked, out of the generated rings, an uncorrelated subset with 
3i topology (positively established with the KNOTFIND algorithm [HI]). Finally a thermodynamic reweighting 
technique [T2] is used to remove the p bias and obtain canonical averages for the observables. 

Figure [l^ illustrates both the range of spherical confinement spanned by the sampled rings as well as the arclength 
covered by the knot. As a first method to identify the knotted portion we look for the arc, r 1; that (i) has the 3i 
topology while the complementary subchain, T2 = r\Ti is unknotted; (ii) cannot be further shortened without losing 
the full-ring 3i topology and finally (iii) can be extended continuously to encompass the whole ring retaining the 3i 
topology and the unknottedness of the complementary arc, T2- The topological state of the open arcs T\ and Y2 is 
established after bridging their termini with the minimally-interfering closure of ref. [17] . 

The above scheme for locating the knot can be seamlessly and robustly applied at all levels of ring compactification. 
As such, it considerably extends the scope of previous investigations where, depending on the free or collapsed nature 
of the rings, the knotted-arc identification had to be made with different and incquivalent methods, including ones 
affecting chain geometry by a rectification procedure or by external constraints such as a shplink [T51 119] . 

Two important facts emerge from Fig.JTJj: for a given ring contour length, L c , the increased confinement is parallelled 
by an increase of the length of the knotted arc, Ik', at the same time, for a given size of the confining sphere one has 
that Ik increases with L c . As a first step towards a comprehensive rationalization of the results we examine the Ik 
versus L c behaviour in the no- and strong-confinement limits, which are shown in Fig. [if. 

In the unconstrained case, it is seen that Ik increases as a power of L c : lk{R c ~> 00) cx L" with a = 0.75 ±0.02. The 
sublinear increase implies that lk/L c vanishes for increasing L Cl resulting in a weak knot localization. This is consistent 
with previous results on knot localization for unconstrained rings both on- and off-lattice where, however, different 
values of a ranging from 0.54 to 0.75 were reported [HI [19]. It also agrees qualitatively with recent experimental 
findings on knotted circular DNA adsorbed on mica surfaces based on single molecule imaging techniques 20 . 

A power-law dependence of Ik on L c provides a good fit of the data for strong confinement too but in this case 
a = 1.09 ± 0.07. This is compatible with a linear dependence of Ik on L c thus implying a full derealization of the 
knot in equilibrated rings subject to severe three-dimensional confinement. To the best of our knowledge this fact has 
not been established before. However, it is worth pointing out that knot derealization has been previously observed 
for 0-collapsed knotted rings |18l 119] . The analogy with our findings is noteworthy since the ring metric properties 
are dictated purely by equilibrium thermodynamics for collapsed rings while spatial constraints are also at play for 
confinement. 

We next turn our attention to the intermediate regime intervening between the no- and strong-confinement limits 
and show that the complex behaviour of Fig. [TJb can be rationalized in terms of a surprisingly simple relationships 
involving the length scales L c , R c and Ik- 
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FIG. 2. (a) Mean squared end-to-end distance, R^ e , of arcs of length s in rings of N = 150 cylinders at different levels 
of confinement. The dashed horizontal line indicates the reference value equal to 6/5i?c while the continuous line shows 
the expected R^ e for an unconstrained, non-self-avoiding Kratky-Porod (KP) chain. The intersection of these reference lines 
determines the nominal deflection arclength, s highlighted by the dropline. (b) Collapse of the same data points as in Fig. [1^ 
obtained by rescaling Ik and L c by s. The slope of the line fitting the data for ^ > 40 is 1.09 ± 0.02. 



To this purpose we first observe that upon increasing confinement the geometric characteristics of the rings resemble 
those of uniformly filled spheres. In particular the eigenvalues of the tensor of gyration of rings become close to the 
eigenvalue of a uniformly-filled sphere (see supporting material). 

The approximately uniform filling of the confining sphere by mildly- or strongly-confined rings is well illustrated 
in Fig. p2a which portrays the mean squared end-to-end separation, i?g e (s) of ring portions of different arclengths, s. 
In fact, the profile of Rl e {s) is noticeably flattened and its plateau value is about 6/5i?^, which is the mean square 
distance of two points inside a spherical volume of radius R c . The value of s at which the plateau value of R\ e is 
reached corresponds to the typical contour length, s, of a discrete worm-like (Kratky-Porod) chain with mean squared 
end-to-end distance equal to 6/5i?^ 



KM = Up 
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s is the nominal arclength of unrestricted ring-portions that are just long enough to hit the sphere boundary and, by 
inverting ([2|, can be expressed in terms of the principal branch of the Lambert W function [21] . 

s/i p = i + y + w(-eM-i-y}) (3) 

where y = 3Rl/5lp. We shall refer to s as the deflection arclength, in analogy with the deflection argument originally 
introduced for polymers in channels [22] and later applied to other contexts, such as crumpled globules [8]. 

The deflection arclength is crucial for rationalizing the results of Fig. [l^ and [if. In particular it is plausible that 
confined rings that experience the same nominal number of deflections at the spherical hull boundary should have 
statistically similar knot lengths once rescaled by s. This conjecture is supported by the scatter plot of Fig. [2}d where 
the rescaled knot length is plotted against the number of nominal deflections L c /s for rings of all considered lengths, 
50 < N < 250, and degrees of confinement. Despite the heterogeneity of the original data sets (see Fig. [!£), the plot 
displays a striking collapse of the data points. A broadening of the curves is seen only in the limit of unconstrained 
rings, where indeed, the deflection length looses meaning. We stress that, because s is calculated deterministically, 
no single adjustable parameter was introduced to obtain the collapse in Fig. (2Jd. The collapsed curve exhibits an 
asymptotic linear trend thus reinforcing the derealization result of Fig. [TJ. 

We conclude the study of the interplay of ring geometry and topology by showing that, by increasing the ring 
compactification, a non-trivial type and degree of entanglement is observed at progressively smaller arclengths. 

To do so we search for the smallest arcs with 3i topology. As illustrated in Fig. [TJ:,d such minimally-long knotted 
portion may differ from the proper knotted arc introduced before because it does not necessarily satisfy criterion (iii) . 
Consequently, it may correspond to an ephemeral knot (because it can be contained in a longer arc with different 
topology, e.g. an unknot) |2"3"j . 

Fig. [3^, illustrates how the length of the shortest knotted arc, l s , behaves for increasing confinement and compares it 
with Ik. It is seen that for no- or midly-constrained rings, when the ring geometrical entanglement is minimal, the two 
measures are in good accord. However, upon increasing confinement they progressively diverge. Notice that l s and 
Ik provide lower and upper limits to the range of arclengths over which non-trivial topological entanglements (after 
arc closure) are observed. Consequently, the divergence of the two metrics indicates that the entanglement cannot be 
described by a single length scale but it displays a multiscale behaviour that amplifies upon increasing confinement. 
Similarly to the case of the proper knot length, Ik, the concept of deflection arclength is useful for rationalizing the 
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FIG. 3. (a) Dependence of knot length, 1^, and the shortest knot length, l s , on confinement for rings of N = 150 cylinders, 
(b) Scatter plot of l s /s versus L c /s for rings of all lengths and degree of confinement. 



interplay of l s , R c and L c . In fact, l s expectedly results from the balance between (a) the increasing probability that 
a knot can be tied with a minimal (nominal) number of deflected segments/cylinders for longer and more confined 
rings and (b) the decreasing probability that the complementary arc is unknotted. 

In the plot of Fig. [SJd it is seen that, for increasing L c , the l s /s data tend to approach a limiting curve that has a 
linear dependence on L c /s. The results indicate that for medium and strong confinement the length of the smallest, 
and possibly ephemeral, knot l s increases approximately linearly both with the L c at fixed R c and with R c at fixed 
L c (more data in supp. mat.). Together with the behaviour of Ik, the results provide a quantitative basis for the 
multiscale character of the ring entanglement that sets in for increasing confinement. 

To summarize, we have used simulations to measure the length of the topologically entangled portion (If.) of knotted 
self-avoiding ring polymers subjected to different degree of spherical confinement. We find that for weak confinement, 
when the geometrical entanglement of the chain is moderate, the knot is weakly localized while for strongly packed 
rings it delocalizes completely. The derealization trend at increasing confinement and different L c follows the neat 
scaling behavior of Fig. [TJ which is rationalised in terms of the deflection arclength. The latter therefore appears to be 
a key quantity to characterize the complex interplay of the geometry and topology in confined polymer rings, which 
reverberates in the above mentioned multiscale entanglement. This observation prompts the question of whether an 
analogous concept and scaling argument can be introduced to characterize and rationalize the entanglement found in 
other systems of densely-packed circular polymers, such as ring melts or theta-collapsed rings [51 [TH1 US] - We wish to 
point out that our results have direct bearings on the relevant problem of viral DNA ejection into infected cells |24) . 
In fact this crucial step of the viral infection is not hindered by the presence, in the packaged genome, of knots because 
the latter are highly delocalised [7J [TT] . 

We thank D. Marenduzzo and A. McGown for useful discussions. We acknowledge support from Italian Ministry 
of Education. 
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SUPPORTING MATERIAL 



A. Size and shape of confined rings 




FIG. 4. The average eigenvalues of the rings gyration matrix, Ax, A2, A3 (largest to smallest) are shown as a function of the 
radius of the confining spherical cavity, R c . Errorbars report the standard deviation around the average value. The inset shows 
the data for high spatial confinement. Going from no- to mild confinement, the average length of the rings principal gyration 
axis is substantially reduced while the other two are much less affected. For example, the initial five-fold reduction of Ai in 
rings of N = 150 cylinders is accompanied by a change of A2 of about 50% and a change of A3 not larger than 10%. The 
latter variation is non-monotonic, indicating a swelling of the mildly confined semi-flexible rings in the plane perpendicular to 
the (squeezed) primary gyration axis. Importantly, the strong shape anisotropy of unconstrained rings is largely lost even for 
mild confinement: already at R c — 5 all average gyration are compatible with the gyration eigenvalue of the uniformly-filled 
confining hull within the estimated errors. 
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B. Multiscale entanglement: detailed comparison of h and l s 
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FIG. 5. Average length of shortest knotted arc (a) for rings of different contour lengths L c and increasing degree of confinement 
1/Rc- Notice that, at a given contour length, l s becomes smaller for increasing confinement. Also, l s increases with L c at fixed 
degree of confinement. The average length of knotted portion, lk, shown in Fig. 1 of the manuscript, is reported in panel (b) 
for ease of comparison. 



(a) 
0.06 



R = 15 



(b) 



R c = 7.5 (c) 



R = 3.5 



1 0.04 



o 0.02 





Is ' 


1 1 


I 
ll 

1 1 




1 1 
1 1 

— 1 1 
1 1 


Ik " 


1 \ 

_ 1 \ 

i ^ 





50 100 150 50 100 150 50 100 150 

lk. Is Ik. Is Ik. Is 



FIG. 6. Probability distribution of the lengths of the shortest knotted arc, l s , and knotted portion, lk, at different levels of 
confinement for rings oi N = 150 cylinders. At no confinement the two distribution are pratically identical, see panel (a). 
However, their difference increases progressively with confinement, see panels (b) and (c). 



C. Comparison with infinitely thin rings 

The robustness and generality of the multiscale character of the entanglement found for self-avoiding rings is 
underscored by the fact that an analogous phenomenology is found for confined infinitely-thin trefoil-knotted rings, as 
reported in the following. The infinitely-thin rings are Kratky-Porod chains: they have the same persistence lengths 
as the self-avoiding ones discussed in the main paper, but have no excluded volume. 

Given the absence of steric hindrance, higher levels of compactifications are achieved compared to self-avoiding 
rings. 



(a) 



Rings with excluded volume 



250 
200 
150 
100 
50 




(b) 



1 


50 


1 

150 - 




■ 75 


175 - 




♦ 100 ► 


200 - 




125 " 


225 _ 






250 . 


— ♦♦♦■♦■♦♦■♦•■♦■'*"*' _ . . ■ • 
. »»»»•»»»■ 


♦ * * 

B a > > 




1.1,1. 


• • • 


♦ ♦ * 

1 



(c) 



250 
200 
150 
100 
50 




250 
200 
150 
100 
50 




0.1 



0.2 



0.3 



0.4 



0.5 



(d) 




250 
200 
150 
100 
50 




50 



100 



150 
N=L c /b 



200 



250 



0.1 



50 



Infinitely thin rings 





1 1 1 . . , 


50 


175 - 






■ 75 ► 


200 




.*,""" 


• 100 . 


225 






' 125 . 


250 






150 












* * J " T ' 


. . « ' A 


♦ ♦ ♦ • 




* * r 


♦ ^ * * _ - ■ ■ 


« ■ ■ 


■ ■ 




■ ■ ■ " " 

• •••••• 




• • 


• 



0.2 



l/Rc 



0.3 



0.4 



0.5 





1 1 1 1 1 1 
R c -^0.5 


I 1 A 






0(= 0.97 +/- 0.02 - 




R c »1 






□ + — 


CX= 0.72+/- 0.02 

1 



100 



150 
N=L c /b 



200 



250 



FIG. 7. (a) and (b): average knot length, as a function of 1/R C for rings with and without excluded volume respectively, 
(c) and (d): Dependence of Ik with L c in the no- and strong- confinement limits for rings with and without excluded volume 
respectively. Note that analogous qualitative behaviour of the two types of models. 
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FIG. 8. Collapse of the data points shown in figure [7b (rings with excluded volume, as in main paper) and[7]l (rings without 
excluded volume) obtained by rescaling Ik and L c by the deflection length, s. The same deflection length as in eqn. 3 of the 
main paper was used in both cases. Note that the slope of the line fitting the data for L c /s > 40 is compatible between the 
two cases. 
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FIG. 9. Top: dependence of knot length, Ik, and the shortest knot length, l s , on confinement for rings of N — 150 segments 
with excluded volume (panel a) and without excluded volume, (panel b). Bottom: scatter plots of l s /s versus L c /s for rings 
with excluded volume (panel c) and without excluded volume (panel d). 



